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Abstract 


In this work, the pulsation analysis is performed on 83 high-amplitude 6 Scuti stars (HADS), which have been 
observed by the Transiting Exoplanet Survey Satellite. The results show that 49 of these HADS show single-mode 
pulsation, 27 of them show radial double-modes pulsation (in which 22 of them pulsate with the fundamental and 
first overtone modes and five of them pulsate with the first and second overtone modes), and seven of them show 
radial triple-modes pulsation (three of which are newly confirmed triple-mode HADS). The histogram of the 
fundamental periods and the ratios between the fundamental and first overtone periods show bimodal structures, 
which might be caused by the stellar evolution in this specific phase. Most of the radial triple-mode HADS have a 
fundamental amplitude of 41—54 mmag, and 50% of them have similar amplitudes of the fundamental and first 
overtone pulsation modes. All these hints require further confirmation not only in observations with more HADS 
samples, but also in theoretical models with suitable treatments of stellar evolution and pulsation. 
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1. Introduction 


ó Scuti stars are a classical type of short-period pulsating 
variable stars, which have periods ranging from 15 minutes to 8 
hr and the spectral classes A-F. In the Hertzsprung-Russell (H- 
R) diagram, they are located on the main sequence (MS), pre- 
MS* or post-MS evolutionary stage at the bottom of the classical 
Cepheid instability strip and are self-excited by the « mechanism 
due to the partial ionization of helium in the out layers 
(Breger 2000; Kallinger et al. 2008; Guenther et al. 2009; 
Handler 2009; Uytterhoeven et al. 2011; Holdsworth et al. 2014; 
Steindl et al. 2022). Most of them pulsate in radial and non-radial 
p-modes (Zong et al. 2015), and some of them also show the 
non-radial g-mode in low-frequency region simultaneously 
(hybrid pulsators, see e.g., Breger & Beichbuchner 1996; 
Bradley et al. 2015; Yang et al. 2021). 

High-amplitude ó Scuti stars (hereafter HADS) are a subclass 
of 6 Scuti stars, which have relatively larger amplitudes 
(AV > 0? 1) and slower rotations (v sini < 30 km s7!) in most 
cases. However, as the accumulation of the HADS samples, these 
classical criteria become unclear now (see e.g., Balona et al. 
2012). Most of HADS pulsate in single or double radial pulsation 
modes (Niu et al. 2013; Alton 2019; Bowman et al. 2021; 
Alton 2022c, 2022; Daszyüska-Daszkiewicz et al. 2022), and 
some of them have three radial pulsation modes (Wils et al. 2008; 


© A detailed guide of asteroseismology for pre-MS stars can be referred to 
Zwintz & Steindl (2022). 


Niu & Xue 2022) or even some non-radial pulsation modes 
(Poretti et al. 2011; Xue & Niu 2020). 

On the aspect of stellar evolution theory, the period variation of 
a single star has significant different values in different 
evolutionary stages. As a result, the observed linear period 
variation rate could be an important criterion to determine the 
evolutionary stage of a star. Based on the times of maximum light 
lasting for decades, some HADS can be considered as normal 
stars evolving into a special evolutionary stage, which can be 
precisely determined according to asteroseismology self-consis- 
tently (see, e.g., Niu et al. 2017; Xue et al. 2018, 2022). These 
results show that HADS should be located in the post-MS 
evolutionary stage. However, some works show that some HADS 
can also be located in the terminal-age main-sequence (TAMS) or 
even MS (see, e.g., Bowman et al. 2021; Sun et al. 2021; Lv et al. 
2022; Yang et al. 2022), although the observed period variations 
are inconsistent with (always significantly greater than) the 
theoretical model predictions in these works. 

These years, a large number of HADS have been monitored 
by the Transiting Exoplanet Survey Satellite (TESS, Ricker 
et al. 2015), whose continuous photometric data provide us to 
study their pulsation properties statistically. 


2. Methods 


We collected 10,731 HADS in the International Variable 
Star Index (VSX), which were then performed cross matching 
with the input catalog of TESS. Finally, we obtained 83 HADS 
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Table 1 
Periods and Amplitudes of the Single-mode HADS 


ID Po(days) Ao(mmag) TESS Sectors 
AD CMi 0.12297 93.4 34* 

AN Lyn 0.09827 53.5 21° 

ASAS J065937-0047.0 0.10882 78.45 33° 

ASAS J015307-5056.5 0.091 87.02 2°, 3, 29, 30 
BL Cam 0.0391 98.76 19° 

BM For 0.05411 74.56 4, 31° 

BO Lyn 0.09335 63.88 21°, 47 

BS Aqr 0.1978 119.28 2°, 29, 42 

CC And 0.1249 47.56 17° 

CY Aqr 0.06104 194.25 42° 

DX Cet 0.10396 56.69 31°, 42, 43 
GP And 0.07868 160.33 17* 

GSC 01951-01755 0.1207 165.6 45, 46° 

GW Uma 0.20317 121.15 21^, 48 

KU Cen 0.08 138.34 37* 

PT Com 0.08211 86.44 22°, 49 

RS Gru 0.14701 154.7 1, 28° 

RY Lep 0.22515 112.61 6, 32°, 33 

SZ Lyn 0.12053 145.8 20°, 47 

TYC 3224-2602-1 0.08866 67.25 16* 

V0337 Ori 0.20129 127.27 43*, 44, 45 
V0358 Mus 0.12177 82.45 11°, 12 
V0367 Cam 0.12161 78.44 19* 

V0398 Uma 0.00134 122.29 15, 21, 22, 41°, 48 
V0411 Sge 0.13669 80.41 40* 

V0451 Dra 0.05581 169.4 14, 20, 21, 40°, 41, 47, 48, 53 
V0467 Dra 0.19768 — 112.31 15, 16, 23, 49, 50 
V0474 Mon 0.13616 60.6 6°, 33 
V0547 Lac 0.09525 96.83 16* 

V0554 Vel 0.10242 93.58 10, 36°, 37 
V0572 Cam 0.08634 117.23 14, 20, 21, 40, 41, 53 
V0575 Lyr 0.14555 57.46 40° 

V0593 Lyr 0.10215 181.1 14, 40°, 53, 54 
V0673 Hya 0.10806 116.7 9", 36 
V0973 Cep 0.08063 118.05 16-18, 19°, 24, 49, 52 
V0974 Oph 0.19102 124 12°, 39 
V1051 Ara 0.11337 . 154.15 12, 39* 
V1307 Sco 0.11703 167.86 12, 39* 
V1338 Cen 0.13093 127.67 11, 38° 
V1421 Cen 0.11417 — 131.93 37* 
V1429 Cen 0.16701 54.69 11? 

V1535 Her 0.09812 62.82 40* 

V2367 Cyg 0.17665 129.02 14?, 15, 40, 41, 54, 55 
V2455 Cyg 0.09421 137.28 15, 16° 
V5505 Sgr 0.0844 136.27 13°, 39 
V6544 Sgr 0.0632 146.7 13°, 27 

XX Cyg 0.13486 227.83 14°, 15-17, 41, 54-57 
YZ Boo 0.10409 119.91 24°, 50, 51 

ZZ Mic 0.06718 123.01 1, 27° 

Note. 


* denotes the Sector used in this work. 


(most of which had 2 minutes cadence flux measurements, 
except some triple-mode HADS) from MAST Portal’ which 
were processed by the TESS Science Processing Operations 
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Center (SPOC; Jenkins et al. 2016). After converting the 
normalized fluxes to magnitudes by utilizing the TESS 
magnitude and removing the long trends in each Sectors, we 
chose the light curves in one Sector (which has the smallest 
value of standard deviation among all the Sectors) of each of 
the 83 HADS to perform the pre-whitening process. 

In the pre-whitening process, the software Period04 (Lenz & 
Breger 2005) was used to perform the Fourier transformations of 
the light curves to search for the significant peaks in the frequency 
spectra from 0 to 150 cd~', until there were no significant peaks 
(S/N>5.6, Zong et al 2018) After removing the alias 
frequencies considering the resolutions and gaps of the data sets, 
we obtained the significant frequencies and their amplitudes for 
each of the HADS. Each of the significant frequencies of an 
HADS was then identified whether it was an independent 
frequency or a harmonics/combination of some independent 
frequencies. Finally, we obtained all the independent frequencies 
(and their amplitudes) for each of the HADS. 

These independent frequencies were identified belonging to the 
radial pulsation modes based on the following relations 
(Stellingwerf 1979): 


0.756 < P,/Py < 0.787, 
0.611 < P5/P, < 0.632, 
0.500 < P3/ Py < 0.525, (1) 


where Po, Pı, P2, and P} are considered as the periods of the 
fundamental, the first overtone, the second overtone, and the third 
overtone pulsation modes respectively. In this work, we would 
strictly follow the relations in Equation (1) to perform the 
identification. Moreover, the independent frequencies which did 
not follow the relations in Equation (1) were considered belonging 
to the non-radial pulsation modes, on which we did not focus in 
this work. For convenience, the fundamental mode, the first 
overtone mode and the second overtone mode are abbreviated as 
F, 10, and 20, respectively." 


3. Results 


In 83 HADS, we find 49 single-mode HADS, 27 radial 
double-mode HADS (22 with F and 10 pulsation and five with 
10 and 20 pulsations), and seven radial triple-mode HADS. 


3.1. Single-mode HADS 


In Table 1, we list the single-mode HADS and their periods and 
amplitudes. In the following statistical analysis, all these single- 
mode HADS are assumed to be pulsating in their fundamental 
modes, which might not be correct in all the cases and needs 
further research based on more information of the stars.” 


3 We did not find quadruple-mode HADS in this work. 


An interesting work (Pietrukowicz et al. 2020) finds that the shapes of the 
light curves are different between the single-mode fundamental and first 
overtone 6 Scuti stars. But this needs a further check and requires clear 
quantitative criteria for application. 
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Table 2 
Periods and Amplitudes of the Radial Double-mode HADS 


ID Po(days) ^ P,(days) —P»(days) P, / Po P3/Pi 
AI Vel 0.111 58 0.08621 0.77264 

BE Lyn 0.095 88 0.07450 0.77705 

FP Cir 0.12675 0.09766 0.77049 

GSC 04 257-00471 0.17379 0.13307 0.76571 

GSC 2583-00504 0.05172 0.03999 0.77321 

GSC 03949-00811? 0.16977 0.13007 0.76615 

NSV 148000? 0.15783 0.12207 0.77342 

NSVS 7293918 0.08854 0.06850 0.77369 

RV Ari 0.09312 0.07195 0.77263 

SX Phe 0.05496 0.04277 0.77821 

V0363 Tra 0.13856 0.10490 0.75710 

V0388 Tel 0.14908 0.11277 0.75649 

0403 Gem? 0.15339 0.11771 ee 0.76735 ee 
V0488 Gem* EE 0.09325 0.07492 ee 0.80342 
V0542 Cam? 0.17477 0.13400 ee 0.76675 ee 
V0703 Sco” 0.14998 0.11521 0.76820 

V0733 Pup 0.22876 0.17424 0.76167 

V0756 Cra 0.10719 0.08216 ee 0.76651 ee 
V0798 Cyg* Ue 0.19478 0.15592 ee 0.80052 
V0899 Car 0.11079 0.08585 ee 0.77486 ee 
V1049 Ara 0.10310 0.08050 0.78079 

V1392 Tau 0.07443 0.05790 ee 0.77789 ee 
V1553 Sco? EE 0.18429 0.14704 ee 0.79783 
V1719 Cyg* Ue 0.26726 0.21374 ee 0.79975 
V2855 Ori 0.05808 0.04483 ee 0.77177 ee 
VX Hya 0.22341 0.17274 ee 0.77322 ee 
VZ Cnc* Ue 0.17837 0.14280 0.80058 


Notes. The ratios of the periods and amplitudes are also listed. 
* denotes the HADS pulsating with 10 and 20. 

5 denotes the HADS whose Ai/Ao 5 1. 

* denotes the Sector used in this work. 


3.2. Double-mode HADS 


In Table 2, we list the 27 radial double-mode HADS and 
their periods and amplitudes, together with the periods and 
amplitudes ratios of the double-modes. What is of interest is 
that there are five HADS pulsating with the 1O and 20 modes 
(V0488 Gem is a newly confirmed such stars) and five HADS 
have larger amplitudes of 1O than F mode, which are denoted 
in Table 2. The origins of these differences should be related to 
the pulsation mode selection mechanisms and worth in-depth 
studies in the future. 


3.3. Triple-mode HADS 


In Table 3, we list the seven radial triple-mode HADS and 
their periods and amplitudes, together with the ratios of the 
periods and amplitudes of different pulsation modes. We also 
note that three HADS have larger amplitudes of 1O than F, 
which might be related to the same cases in double- 
mode HADS. 

Up to now, only 15 HADS are confirmed to pulsate in the 
first three radial modes (the fundamental mode plus the first and 


Ao(mmag)  A,(mmag) A (mmag)  Aj/Ao  Ao/Ai TESS Sectors 
93.61 65.79 0.70 : 34°, 35 
114.65 2.59 0.02 Ut 21° 
78.34 2.84 0.04 12, 38° 
102.24 52.71 0.52 16, 17, 18, 24* 
54.61 18.31 0.34 245, 51, 52 
6.21 8.27 1.33 16, 17, 41°, 55 
54.84 76.57 1.40 1, 2°, 28 
118.69 8.54 0.07 20°,44-47 
135.90 39.67 0.29 42, 43° 
138.18 31.57 0.23 2°, 29 
140.34 3.06 0.02 12, 39* 
56.16 36.95 0.66 13°, 27 
42.42 43.94 Ur 1.04 43°, 44, 45 

Zn 129,53 0.80 Ur 0.01 33* 
71.09 84.71 1.19 Uy 19* 
49.59 73.41 1.48 Ut 39* 
107.69 47.62 0.44 U 7, 8, 34° 
104.07 30.02 Us 0.29 Ut 13* 

Us 144.60 15.60 Ut 0.11 14, 40, 41°, 54 
99.22 31.13 Us 0.31 U 10, 11, 37* 
118.82 8.97 0.08 Zn 39* 
81.86 43.89 e 0.54 Us 5* 

EE 52.45 25.87 Uer 0.49 12, 39* 

Us 96.27 10.31 Ut 0.11 155, 16, 55 
88.59 24.70 Us 0.28 U 65, 33 
92.00 74.71 Zn 0.81 Us 8, 35° 

Us 125.09 50.04 Ut 0.40 7, 34, 44, 45, 46° 


second overtone) in the galaxy: V829 Aql (Handler et al. 
1998), GSC 762-110 (DO CMi) (Wils et al. 2008), GSC 
03 144-595 (Mow et al. 2016), GSC 08 928-01300 (Yang et al. 
2021) V761 Peg (Kazarovets et al. 2020), two cases in 
Khruslov (2014) (V0803 Aur and V1647 Sco), four cases in the 
OGLE project (OGLE-GD-DSCT-0021, OGLE-GD-DSCT- 
0033, OGLE-GD-DSCT-0048, and OGLE-GD-DSCT-0049), 
and four cases in Khruslov (2022) (GSC2.3 NBY9001324, 
GSC2.3 NAVT000282, GSC2.3 S5SR001526, and GSC2.3 
S4NMO025687). In addition, two triple-mode HADS are 
confirmed in the Large Magellanic Cloud (OGLE-LMC- 
DSCT-0927 and OGLE-LMC-DSCT-2345) (Poleski et al. 
2010). These stars exhibit different pulsating properties 
compared with other HADS, which would help us to under- 
stand the selection mechanism of pulsation modes deeply (Niu 
& Xue 2022). 

In this work, three radial triple-mode HADS are newly 
confirmed: ASAS J094303-1707.3, V1384 Tau and V1393 
Cen, and their light curves and frequency spectra are shown in 
Figures 1, 2, and 3, respectively. The continuous photometric 


Table 3 

Periods and Amplitudes of the Radial Triple-mode HADS 
ID Po(days) P(days) P2(days) P,/Po P2/Po Ao(mmag) A1(mmag) A»(mmag) A,/Ao A2/Ao TESS Sectors 
ASAS J094303-1707.3* 0.09918 0.07652 0.06134 0.77149 0.61852 53.54883 61.60368 2.01469 1.15042 0.03762 8,35° 
GSC 762-110* 0.19450 0.14862 0.11909 0.76415 0.61228 48.988 91 49.01288 25.00426 1.00049 0.51041 71800 s), 33*(600 s) 
GSC 03144-595 0.20367 0.15548 0.12445 0.76338 0.61104 50.66162 49.35125 11.25300 0.97413 0.22212 14(1800 s), 15*(1800 s) 
GSC 08928-01300^ 0.07573 0.05881 0.04750 0.77656 0.62720 76.66062 0.55071 0.91022 0.00718 0.01187 4,5,7-11,27,28,31,34,35,37,38° 
V0803 Aur 0.07106 0.05503 0.04439 0.77448 0.62473 49.48538 20.77060 4.49043 0.41973 0.09074 43°,44,45 
V1384 Tau 0.13980 0.10739 0.08589 0.76820 0.61439 52.75469 49.93025 4.40794 0.94646 0.08356 5° 
V1393 Cen* 0.11778 0.09083 0.07285 0.77118 0.61850 41.30245 70.78809 5.65266 1.71390 0.13686 11,38° 


Notes. The ratios of the periods and amplitudes are also listed. 

è denotes the HADS whose A;/Ao > 1. 

> denotes the HADS is a 6 Scuti and Dor hybrid star. 

* denotes the Sector used in this work; GSC 762-110 and GSC 03144-595 do not have the data of 120s exposure, we used the 600 and 1800 s data instead. 
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Figure 1. Light curves and frequency spectrum of ASAS J094303-1707.3 observed by TESS. 
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Figure 2. Light curves and frequency spectrum of V1384 Tau observed by TESS. 
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data provide us an opportunity to study the nature of them in 
the near future. 


4. Discussions and Conclusions 


Although the number of the samples is not that large, the 
statistical analysis of the pulsation properties of all the 83 
HADS in Tables 1, 2, and 3 could provide us some hints about 
the different origins of the three types. 

Figure 4 shows the distribution of Po (period of the 
fundamental pulsation mode) of the 78 HADS (except the 
radial double-mode HADS pulsating with 1O and 20), in 
which we can see that the most values of Po are concentrated in 
the range of 0.08 to 0.13 days. In detail, there is a small peak at 
about 0.2 days, which might be related to the two possible 
evolutionary stages of HADS (MS and post-MS). 

Figure 5 shows the distribution of A; (amplitude of the 
fundamental pulsation mode) of the 78 HADS (except the radial 
double-mode HADS pulsating with 10 and 20), in which we 
can see that the most values of A, are concentrated in the range 
of 30-130mmag. The samples decrease slowly when 
Ag 130mmag while drops sharply when Ap < 30 mmag 
(which might be caused by the selection criterion of the HADS 
samples). All the triple-mode HADS have the Ag in the range of 
30-80 mmag, if we consider the particularity of GSC 08 928- 
01300 (a 6 Scuti and y Dor hybrid star, Yang et al. 2021), 
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Figure 3. Light curves and frequency spectrum of V1393 Cen observed by TESS. 


all the other triple-mode HADS have the values of Ap from about 
4] to 54 mmag. The low-amplitude of these HADS can be 
explained by the fact that more pulsation modes have to share the 
driving energy in the ionization zones, while the concentrated 
distribution of them might be a feature of this type of stars rather 
than a coincidence, which needs more samples to confirm. 

Figure 6 shows the distribution of P, /Po (periods ratio of the 
first overtone by the fundamental pulsation mode) of the 29 
HADS (except the single-mode HADS and double-mode 
HADS pulsating with 10 and 20), in which that of the 
double-mode shows an obvious bimodal structure, while that of 
the triple-mode shows a relative smooth distribution from about 
0.763 to 0.778. 

Figure 7 shows the distribution of A, /Ag (amplitudes ratio of 
the first overtone by the fundamental pulsation mode) of the 29 
HADS (except the single-mode HADS and double-mode 
HADS pulsating with 10 and 20), in which we can see that 
the most values of A,/Ag of double-mode HADS are 
concentrated in the range of about 0 to 0.5. What is interesting 
is that, over 6046 of the triple-mode HADS have the values of 
A1/Ao around 1.0 (from about 0.8 to 1.2), which could be 
ascribed to the coupling between the fundamental and first 
overtone modes of these stars. Such as the case in the triple- 
mode HADS KIC 6382916 (GSC 03144-595), the three 
pulsation modes are related by the resonance relationship 
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Figure 5. Distribution of Ap of the 78 HADS. 


2f, + Aw = fo + fa (where Aw is the rotation frequency of the 
star), which couples the fundamental and first overtone modes 
together (Niu & Xue 2022). 

In Figure 8, we plot the Petersen diagram of the 22 HADS 
(which have fundamental and first overtone pulsation modes 
except the single-mode HADS and double-mode HADS 
pulsating with 1O and 20), together with the evolutionary 
tracks of single stars calculated by MESA (Paxton et al. 
2011, 2013, 2015, 2018, 2019) and GYRE (Townsend & 


Teitler 2013; Townsend et al. 2018; Goldstein & Town- 
send 2020), with the solar element abundances: X — 0.7438, 
Y — 0.2423, and Z — 0.0139 (Asplund et al. 2021); the mixing- 
length parameter: ayyr=1.89 (Niu et al. 2017); the 
exponential scheme overshooting (Herwig 2000) whose 
parameter depends on the stellar mass M (in solar masses): 
fov = (0.13M — 0.098)/9.0 (Magic et al. 2010) from MS to 
post-MS with masses from 1.5 to 2.5 Mo with a step of 0.1 Mo. 
As a comparison, the evolutionary tracks are also plotted in the 
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Figure 6. Distribution of P,/Po of the 29 HADS. 
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Figure 7. Distribution of A,/Ao of the 29 HADS. 


H-R diagram. The period change rate of the fundamental 
mode (which can be treated as a marker of the stellar 
evolutionary rate of the star) is also plotted in different 
regions: 10 5 yr ! < | Py /Pol<1077 yr-!, 1077 yr“! < |É /Pp 
|< 10$ yr^!, and [Py /Po| > 107 yr. 

In each of the evolutionary tracks, the round evolutionary 
phase corresponds to the overall contraction phase after MS 
(i.e., the commonly known “hook” in the H-R diagram), 


which is a quite rapid evolutionary phase in the range from 
MS to post-MS (Aerts et al. 2010). What is interesting is that 
these most rapidly evolving states in each of the evolutionary 
tracks start when P,/Po ~ 0.77, which is clearly represented 
in the Petersen diagram of Figure 8. As a result, if we assume 
that the HADS are normal stars evolving in the MS to post- 
MS phase, the distribution of their Pı /Po would have a gap at 
about 0.77 and show a bimodal structure, which is consistent 
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Figure 8. Petersen diagram (upper panel) of the 29 HADS and the relevant H-R diagram (lower panel) of the evolutionary tracks. The blue dots represent the double- 
mode HADS, and the red dots represent the triple-mode HADS. The colored thin lines represent the evolutionary tracks from MS to post-MS with masses from 1.5 to 
2.5 Ms with a step of 0.1 Mo, while the thick gray, orange, and red lines represent the evolutionary states whose |P)/P| fall into the range 
10-8 yr-! < |Py/Pol < 1077 yr-!, 1077 yr-! < [Py /Po| < 1076 yr-!, and |B) /Po| > 1079 yr-!. 


to the result shown in Figure 6.'° Tt gives us a hint that the 
pulsation properties of HADS might be related to their 


10 Although the evolutionary tracks are simply constructed based on the solar 
element abundances, the tracks with different element abundances have similar 
round evolutionary phases (even though they will have overall panning). If we 
combine all these tracks together, there still exists a most rapidly evolving state 
in the Petersen diagram around P; /Po ~ 0.77. 


evolutionary phases, which should be tested and confirmed 
with more samples in the future. 
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